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We empirically determine effective temperatures and bolometric luminosities for a 
large sample of nearby M dwarfs, for which high accuracy optical and infrared photom- 
etry is available. We introduce a new technique which exploits the flux ratio in different 
bands as a proxy of both effective temperature and metallicity. Our temperature scale 
for late type dwarfs extends well below 3000 K (almost to the brown dwarf limit) and 
is supported by interferometric angular diameter measurements above 3000 K. Our 
metallicities are in excellent agreement (usually within 0.2 dex) with recent determi- 
nations via independent techniques. A subsample of cool M dwarfs with metallicity 
estimates based on hotter Hipparcos common proper-motion companions indicates 
our metallicities are also reliable below 3000 K, a temperature range unexplored un- 
til now. The high quality of our data allow us to identify a striking feature in the 
bolometric luminosity versus temperature plane, around the transition from K to M 
dwarfs. We have compared our sample of stars with theoretical models and conclude 
that this transition is due to an increase in the radii of the M dwarfs, a feature which 
is not reproduced by theoretical models. 
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1 INTRODUCTION 

Low-mass dwarfs are the dominant stellar component of the 
Galaxy and have been employed in a variety of Galactic 
studies: tracing Galactic disk kinematics (e.g. Hawley, Gizis 
& Reid 1996; Gizis et al. 2002; Bochanski et al. 2005, 2007), 
studying the stellar age-velocity relations (West et al. 2006), 
investigating Galactic structure (e.g. Reid et al. 1997, Ker- 
ber et al. 2001, Pirzkal et al. 2005), and the Galaxy's mass 
and luminosity (e.g. Hawkins & Bessell 1988; Kirkpatrick 
et al. 1994; Gould, Bahcall & Flynn 1997, 1998; Zheng et 
al. 2001, 2004). An increasing number of M dwarfs are now 
known to host exoplanets (e.g. Delfosse et al. 1998; Butler et 
al. 2004; Bonfils et al. 2007, Udry et al. 2007). The determi- 
nation of accurate fundamental parameters for M dwarfs has 
therefore relevant implications for both stellar and Galactic 
astronomy. 

Observationally, the spectra of these stars are marked 
by the presence of strong molecular absorption features, in 
either optical (e.g. TiO and VO) and infrared regions (e.g. 
H2O and CO). Molecular lines blend with all other lines and 
create a pseudo-continuum, rendering all spectral analysis 
difficult (e.g. Gustafsson 1989) . However, recent advances in 
model atmospheres of low-mass dwarf stars (Hauschildt et 
al. 2003, Brott & Hauschildt 2005), have boosted the num- 
ber of studies deriving accurate metallicities for M dwarfs 



(Woolf & Wallerstein 2005, 2006; Bean et al. 2006a, 2006b). 
Modeling the internal structure, atmospheric properties and 
magnetic activity of M dwarfs (e.g. Burrows et al. 1993, Al- 
lard et al. 1997, Baraffe et al. 1998, Hauschildt, Allard & 
Baron 1999, Alfred et al. 2006, Reiners & Basri 2007) also 
presents ongoing theoretical challenges. For a small num- 
ber of nearby M dwarfs, interferometry is currently provid- 
ing direct angular diameter measurements (Lane, Boden & 
Kulkarni 2001; Segransan et al. 2003; Berger et al. 2006) 
which confirm a large discrepancy between the predicted 
and observed radii, as has been noted in eclipsing binaries 
with M-type components (see Ribas 2006 for a review). 

In this paper, we empirically determine the effective 
temperatures and the bolometric luminosities for more than 
340 M dwarfs. This work is an extension of our previous 
study on G and K dwarfs, to which we applied the InfraRed 
Flux Method (IRFM, Casagrande, Portinari & Flynn 2006). 
The effective temperature and the bolometric luminosity 
scales we derive are accurate at a level of a few percent and 
are supported by interferometric angular diameters above 
~ 3000 K. 

We find in this study that below about 4000 K the 
monochromatic to bolometric flux ratio in different bands 
is a proxy of both effective temperature and metallicity. By 
exploiting this feature, we are able to derive not only T e s, 
but also the metallicities of our M dwarfs, which are found 
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to be in very good agreement (usually within 0.2 — 0.3 dex) 
with those inferred using the Bonfils et al. (2005) calibra- 
tion or directly measured from Woolf & Wallerstein (2005, 
2006). The technique we propose also looks promising for 
stars much cooler than those explored in the aforementioned 
studies. 

We find considerable structure in the temperature- 
luminosity plane, especially around the transition between 
K and M dwarfs. Our study circumstantially confirms pre- 
vious works which indicate that the radii of M dwarfs are 
larger by about 15% than model predictions. We also find 
strong evidence that this discrepancy, clearly observed in M 
dwarfs in eclipsing binary systems, is also present in nearby 
disk M dwarfs. 

The paper is organized as follows. In Section [2] we de- 
scribe our sample of M dwarfs and in Section[3]we compare it 
with the Phoenix model atmosphere in the two-colour plane. 
We then review the IRFM and present our new technique 
for estimating effective temperatures and metallicities be- 
low 4000 K in Section [4] Our proposed metallicity, bolomet- 
ric luminosity and effective temperature scales along with 
the comparison with other empirical determinations are dis- 
cussed in Section [5] [6] and [7] respectively. In Section [8] we 
analyze the stars with good Hipparcos parallaxes in the HR 
diagram. We find a strong feature which marks the transi- 
tion from K to M dwarfs and which is due to an increase 
in the observed stellar radii. We briefly discuss possible rea- 
sons, including the effect of magnetic fields and molecular 
opacity. We finally conclude in Section [9] 



2 THE SAMPLE 

Our basic sample consists of 343 nearby M dwarfs with high 
quality optical and infrared photometry. We describe the 
sample selection in this Section. 

2.1 Johnson— Cousins photometry 

In recent years major efforts have been devoted to obtaining 
high accuracy Johnson-Cousins photometry for cool stars 
(Kilkenny et al. 1998, 2007; Koen et al. 2002). One third 
of our sample is built from the extensive work of Koen 
et al. (2002) who presented homogeneous and standardized 
UBV(RI)c photometry for more than 500 M stars, half of 
which are main sequence dwarfs. Variability has a very high 
incidence among cool stars; however the existence of ~ 100 
or more Hipparcos measurements per star spread over sev- 
eral years, together with the excellent temporal stability of 
the magnitude scale, permits the detection of variability at 
the level of few hundredths of a magnitude. The red stan- 
dards provided by Koen et al. (2002) are all non-variables 
in this sense (i.e. the Hipparcos variability flag is not set); 
Koen et al. (2002) also provide the SIMBAD and CCDM 
(Dommanget & Nys 1994) classification for variability and 
double/multiple stars. We have discarded all star having 
those labels, for a total final sample 128 M dwarfs. Besides 
Johnson-Cousins photometry accurate to 0.01 mag or bet- 
ter, all these stars have Hipparcos parallaxes better than 
15% and are all closer than 60 pc. 

Another accurate source of Johnson-Cousins photom- 
etry for M dwarfs is Reid et al. (2002, 2003, 2004). Alto- 



gether, they provided new V(RI)c (and B for the bright- 
est sources) measurements for more than 370 NLTT stars. 
We have used the SIMBAD classification plus the essential 
notes as given in the aforementioned papers to remove dou- 
ble, variables, flares, possible misclassification and stars with 
nearby companions (or background stars) that could affect 
the photometric measurements. The photometry is accurate 
and consistent with the standard photometric system to bet- 
ter than 1% (Reid et al. 2002) and with a typical uncertainty 
less than 2% even for single night observations (Reid et al. 
2003, 2004). Furthermore, all observations are done with the 
same instrument and reduced using identical methods, sim- 
ilar to those described in Kilkenny et al. (1998) and Koen et 
al. (2002). Reid et al. (2002, 2003, 2004) also provide photo- 
metric distances for all the stars: although such estimates do 
not provide definitive distance measurements, we have used 
them to ensure that all the stars chosen are closer than 50 
pc and their absolute magnitudes are consistent with those 
expected from dwarf stars. Altogether, we have retained 94 
stars with BV(RI)c and 157 with V(RI)c photometry. 

We also took a few very red dwarfs from Henry et al. 
(2004) who measured V(RI)c colours using standards from 
Graham (1982), Bessell (1990a) and Landolt (1992). These 
stars are also generally fainter and the accuracy is somewhat 
lower, with a typical uncertainty of ±0.03 mag in each band 
(Henry et al. 2004). We have selected 11 such stars, all closer 
than 40 pc. Finally, another source of V(RI)c photometry 
for very red stars is from Bessell (1991), from which we took 
5 stars. For the stars in Henry et al. (2004) and Bessell 
(1991) we have used the SIMBAD classification to avoid 
flares or variables, although in this case the most affected 
bands would be the blue ones (U and B) which we do not 
have for these stars. 



2.2 Near— Infrared photometry 

All the stars presented in Section [27T] have JHKs photome- 
try from 2MASS. In what follows we will use only stars with 
total photometric errors (as given from 2MASS) smaller 
than 0.10 mag (i.e. "j_" + "h_" + "k_msigcom"< 0.10), thus 
reducing our final sample to 343 stars. The typical errors in 
J and Ks are around 0.023 mag, whereas a slightly larger 
uncertainty (0.032 mag) affects H band photometry. 

For most of the Koen et al. (2002) stars, excellent SAAO 
JHK photometry is also available (Kilkenny et al. 2007) : 
this additional photometry is very valuable to check the ac- 
curacy of the 2MASS zero-points and the dependence of the 
proposed temperature scale on the adopted absolute calibra- 
tion. As we show in Appendix B, the SAAO IR photometry 
confirms the adopted 2MASS zero-points and absolute cali- 
bration. 



2.3 Metallicities 

Measuring metallicities for M dwarfs is still challenging. 
With decreasing temperature, the spectra show increasingly 
abundant diatomic and triatomic molecules. The molecu- 
lar bands complicate the calculation of stellar model atmo- 
spheres and cause line blends, making it difficult to estimate 
the true continuum and to measure atomic line-strengths 
over large regions of the visible spectra. 
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Figure 1. Synthetic optical and IR colours from the Phoenix models compared to our 343 sample stars of Section [2] Lines correspond 
to models with [M/H] equal +0.5 (black), +0.0 (red), —0.5 (cyan), —1.0 (green) and —1.5 (blue). Since metallicities for M dwarfs are 
usually not available or are very uncertain, we have not used any metallicity bin for the stars (open squares). Typical error bars are 
shown only for certain indices; in the other bands error bars are comparable or smaller than the size of the plotting symbols, with very 
red stars (from Henry et al. 2004) which are somewhat less accurate. 



Major advances in the field have been recently obtained 
using high-resolution spectra to measure equivalent widths 
of atomic lines in regions not dominated by molecular bands 
(Woolf & Wallerstein 2005, 2006). Bonfils et al. (2005) have 
measured the metallicity in 20 binary systems, having an M 
dwarf secondary and a warmer primary (for which a metal- 
licity is much more readily obtained). They combined their 
results with the abundances measured from Woolf & Waller- 
stein (2005) to calibrate an absolute K$ band luminosity 
versus colour relation (Mk s and (V — Ks)) as a function 
of metallicity. This results in a metallicity relation for M 
dwarfs, but since it depends on absolute magnitude, we ap- 
ply it only to those of our stars with accurate parallaxes 
available from Hipparcos. Altogether in our sample there 
are 118 stars with accurate 2MASS photometry (see Section 
I2.2|l . Hipparcos parallaxes better then 15% and within the 
range of applicability of the Bonfils et al. (2005) calibration. 
The formal accuracy of the relation is ±0.2 dex, however the 
uncertainty in parallaxes introduces an additional error that 
in the worst case can be of the same magnitude. Therefore, 
we expect these metallicities to be accurate to 0.2 — 0.3 dex. 

In Section [4] we will use these metallicities to calibrate 
our method to estimate [M/H] for the rest of the M dwarfa3- 

1 The calibration of Bonfils et al. (2005) returns [Fe/H], whereas 
the model atmospheres are given for the total heavy elements 



3 THE PHOENIX MODEL ATMOSPHERES 

We apply an extension of the IRFM we developed for G 
and K dwarfs (Casagrande et al. 2006) to M dwarfs. While 
most of the bolometric flux of the stars is emitted in the 
optical and infrared and is covered by our observations, as 
in the earlier work, we use model atmospheres to estimate 
the small part (typically less than 20%) of the flux which is 
outside our observational bands. 

Pioneering work on M dwarfs model atmospheres trace 
back to Tsuji (1969) and Auman (1969). The inclusion of 
sophisticated physics became available with the work of 
Mould (1975, 1976) and has steadily continued with Allard 
(1990), Kui (1991), Brett & Plez (1993), Allard & Hauschildt 
(1995), Brett (1995a, 1995b), Tsuji, Ohnaka & Aoki(1996), 
Hauschildt et al. (1999). 

We use the most recent grid of model atmosphere pub- 
licly available at the Phoenix project's websit^S The models 
cover a range of parameters far wider than that needed for 
the present work: 2000 < T cft < 10000 K, -0.5 < log( 5 ) sC 

fraction [M/H]. For low values of alpha— enhancement, the differ- 
ence between the two is negligible, particularly since metallicity 
measurements in M dwarfs are still uncertain. In the rest of the 
paper we will refer to [M/H], although this is in practice [Fe/H] 
when we refer to empirical measurements. 
2 ftp. hs. uni-hamburg. de /pub/ outgoing/phoenix / GAIA / 
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5.5 and -4.0 < [M/H] < +0.5. Below 7000 K the grid is 
given in steps of 100 K in effective temperature and 0.5 dex 
in metallicity. The molecular line lists include about 700 mil- 
lion molecular lines, 15-300 million of which are typically 
selected in a model. The equation of state is an extension 
of that used in Allard & Hauschildt (1995). For the coolest 
models the dust is assumed to form and to immediately 
rain out completely below the photosphere ( "cond" models) 
so that it does not contribute to the opacity. Full details 
are available in Brott & Hauschildt (2005) and references 
therein. 

Since we are working with dwarf stars, we assume 
Iog(<?) = 5.0 throughout. This assumption is in agreement 
with the values of \og(g) determined from other techniques 
(Segransan et al. 2003; Berger al. 2006; Bean et al. 2006a). 
As we will see later, a change of ±0.5 dex in the assumed 
surface gravity implies only minor differences in the derived 
parameters. 

The Phoenix models also include variations of a- 
elements for each metallicity. We have chosen to use models 
with no a-enhancement; in any case the use of a-enhanced 
models in dwarfs of earlier spectral type does not change the 
results significantly (Casagrande et al. 2006). There are indi- 
cations that M dwarfs follow the same [a/Fe] vs. metallicity 
as measured in FGK dwarfs (Woolf & Wallerstein 2005) and 
since our sample is limited to the solar neighborhood we do 
not expect significant signs of a-enhancement. 



3.1 Colour— Colour plots 

Testing of synthetic model atmospheres is normally done 
by comparing observed and modelled spectral energy distri- 
butions for a range of wavelengths and spectral types (e.g. 
Tinney, Mould & Reid 1993; Brett 1995b; Leggett et al. 
1996, 2000, 2001; Burgasser, Cruz & Kirkpatrick 2007). Here 
we simply compare synthetic and empirical colours, showing 
that with good accuracy data right across the optical and in- 
frared, photometry provides an excellent tool to test model 
atmospheres (e.g. Bessell, Castelli & Plez 1998). The inter- 
ested reader can find e.g. in Leggett et al. (2000, 2001) and 
Burgasser et al. (2007) a thorough discussion of the compar- 
ison and pitfalls between observed and modelled spectra for 
M dwarfs and subdwarfs. 

Synthetic colours are computed from the model at- 
mospheres using prescriptions very similar to those in 
Casagrande et al. (2006), as we discuss in more details in 
Appendix B. In Figure [T] we plot various colour indices as 
function of V — Ks ■ It is clear from the figure that the mod- 
els have problems in U — B, being offset with respect to the 
observations by 0.1 — 0.2 mag. Such an offset is unlikely to be 
entirely due to uncertainties in standardize the U band filter 
(see Appendix B) and/or the zero-points (as these amount 
to be at most 0.01 — 0.02 mag) and most likely reflects in- 
adequacies in modeling the UV spectral region in cool stars. 
The models behave considerably better in the other optical 
and infrared bands, although J — H and J — Ks appear off- 
set when going to dwarfs redder then V — Ks > 5.5 which 
corresponds to T E s < 3000 K (compare with Figure [9}. 
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Figure 3. Same as Figure [T] but as function of other colour indices. 



Alternatively to V — Ks, another excellent temperature 
indicator in cool stars is Ic — J- In fact, Figure [2] looks very 
similar to Figure [1] but the stars are now distributed over a 
shorter baseline. 

In the left panel of Figure [3] we show infrared indices as 
function of (7? — I)c, which is sensitive to decreasing spec- 
tral type for M dwarfs (Bessell 1991). The models are in fair 
agreement with the data, but get worse going to the reddest 
(R — T)c. Although at the very red end observed colours 
might be slightly less accurate (Section 12. \\ , Bessell (1991) 
has shown that the large spread in observed colours for the 
latest M dwarfs is real. The reason is likely to be that at 
such cool temperatures M dwarfs show the effects of dust at 
bluer wavelengths and have slightly different spectra with 
stronger hydride bands and weaker TiO and VO bands. In 
the purely infrared colour planes of Figure UJ we are work- 
ing in narrow colour ranges and observational errors become 
prominent in the comparison. To help establish the trends 
we have also checked how these planes look when our sample 
of M dwarfs is complemented with the earlier one for G and 
K dwarfs from Casagrande et al. (2006). At bluer colours the 
data show a turnover in J — H and J — Ks as function of 
decreasing temperature (i.e. increasing H — Ks and Ic — J)- 
The models partly predicted this feature, which is the result 
of the sensitivity of infrared colours to the photospheric gas 
pressure (Mould 1976) as well as to the occurrence of H2O 
bands. At very low effective temperatures, the data suggest 
a flattening and a possibly a rising in J — H and J — Ks, 
whereas the models decrease steadily. The rising of the ob- 
served colours in the H — Ks vs. J — H plane is confirmed by 



similar plot using dwarfs much cooler than we have here (e.g. 
Reid et al. 2001; Burgasser et al. 2007). This mismatch be- 
tween data and models was already noticed in other models 
by Brett (1995b), and essentially means that in the models 
H and Ks magnitudes become progressively fainter than J 
magnitude as the effective temperature decreases. 

Overall, the Phoenix models show fair agreement with 
the data in various bandpasses, although inadequacies still 
persist, especially at the coolest temperatures where the 
dust needs to be properly incorporated. The synthetic 
colours also show a large spread with metallicity for decreas- 
ing T ff: since our sample is limited to the solar neighbour- 
hood, we expect our M dwarfs share a distribution similar to 
that observed in GK dwarfs, i.e. with most of the metallici- 
ties between —0.5 and 0.0. To this respect, the coolest mod- 
els are considerably offset with respect to the position of the 
stars which are encompassed by super-solar and solar lines, 
rather than solar and sub-solar as one would expect from the 
argument mentioned above. The serious discrepancies in U 
band have little impact on the present work of calibrating 
the bolometric luminosities of the stars, since so little flux 
is emitted in the U band. For the reddest stars, which cor- 
respond to effective temperatures below 3000 — 2900 K (see 
Figure [9} the models do show some problems, especially in 
the infrared. Since our technique partly relies on model at- 
mospheres, this means that at the very cool end the results 
we present in Section [5] to [7] are still open to refinements. 
Looking at Figures [T] to [3] is obvious that by simply compar- 
ing observed and theoretical colours one would deduce dif- 
ferent stellar parameters (effective temperatures and metal- 
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licities) depending on the colour index adopted for the com- 
parison. The technique we present in the next Section is less 
affected by such inconsistencies. In fact, we will use model 
atmosphere only to estimate the flux outside our multi-band 
photometry, i.e. the flux in the blue and red tails of the spec- 
tra. This estimate should be rather accurate as long as model 
atmosphere reproduce the overall spectral energy distribu- 
tion, even though specific bands might have problems. No- 
tice though, where theoretical models do show limitations, 
like in the optical and near-infrared, we use anyway the ob- 
served colours. In addition, for estimating the metallicities 
we will calibrate our technique on other empirical measure- 
ments (see Section [5} , and this should keep under control 
the deficiencies that still affect theoretical models. Nonethe- 
less, it would be too optimistic to believe we are not affected 
by the inaccuracies in the models. Future improvements on 
the theoretical side will certainly benefit to our technique, 
especially below 3000 K. 



4 THE MULTIPLE OPTICAL-INFRARED 
TECHNIQUE 

In our previous paper we used multi-band photometry to 
implement the IRFM and we derived effective temperatures, 
bolometric luminosities and angular diameters for a set of G 
and K dwarfs (Casagrande et al. 2006). It is therefore nat- 
ural to ask whether the same technique can be successfully 
applied to M dwarfs. 

Although the underlying idea of the IRFM is still 
valid when going to effective temperatures cooler than ~ 
4000 K, some caveats exist. Here, we generalize and extend 
our temperature scale to dwarf stars much cooler than in 
Casagrande et al. (2006) to which the reader can refer for 
an introduction to the formalism and details on the com- 
putation of the bolometric and monochromatic flux from 
multi-band photometry. The extension of the method pre- 
sented here concerns the computation of T e g, which is now 
done by using the fluxes in both optical and infrared bands 
(Section • For this reason we call our method the Multi- 
ple Optical-Infrared Technique (MOITE). As in Casagrande 
et al. (2006), the effective temperatures we derive depend on 
very few basic assumptions, namely the adopted Vega abso- 
lute calibration and zero-points (see also Appendix B). The 
dependence on the adopted grid of model atmospheres is 
also not so crucial since most of the bolometric flux (usually 
around 80 percent) is recovered from our multi-band pho- 
tometry. The MOITE proves to be also sensitive to metal 
content in M dwarfs, as we discuss in Section [4.31 



4.1 The IRFM in brief 

The basic idea of the IRFM (Blackwell & Shallis 1977; 
Blackwell, Shallis & Selby 1979; Blackwell, Petford & Shal- 
lis 1980) is to compare the ratio between the bolomet- 
ric flux TboI (Earth) and the infrared monochromatic flux 
T\ (Earth), both measured at Earth (the so called observa- 
tional R i, a factor) to the ratio between the surface bolomet- 
ric flux (crTeff 4 ) and the surface infrared monochromatic flux 
J~\ (model) , predicted from model atmospheres. The ratio of 
the last two quantities defines the theoretical Rtheo factor. 
From this ratio T e g can be computed iteratively as follow: 



Teff.n = 



J^a (model) (n _ 1) J' Soi (Earth) (n _ 1} 
ct.Fa (Earth), 1} 



(1) 



where the effective temperature determined at the n th 
iteration depends on the effective temperature determined 
at the (n — l) th iteration and which is used to improve the 
estimate of the quantities on the right hand side of equation 
fTJ. In the IRFM, more than one infrared band is usually 
used (i.e. Rtheo is computed for each band), and the proce- 
dure described here is applied to each band separately. At 
each iteration the average T e g tXl obtained from all the in- 
frared bands is then computed and the procedure is iterated 
until the effective temperature converges to a final value. 

The IRFM has been traditionally applied to stars of 
K or earlier spectral type, to derive effective temperatures 
approximately above 4000 K. Qualitatively, above this tem- 
perature, spectra roughly behave like black-body curves in 
the infrared, so that in this region the spectra can be de- 
scribed by the Rayleigh-Jeans law and the ratio between the 
bolometric and monochromatic flux depends on the effective 
temperature to some power, with little or no metallicity de- 
pendence (Figure 3}. 



4.2 From the IRFM to the MOITE: a top level 
description of the technique 

We now generalize the technique presented in the previous 
Section to effective temperatures cooler than 4000 K using 
both optical and infrared bands with the MOITE. We aim 
to give here a qualitative description of the underlying idea, 
leaving the full technical details to Appendix A. 

Below approximately 4000 K, molecular absorption and 
flux redistribution become very important and significantly 
change the original continuum shape, making any type of 
qualitative black-body description to an M dwarf spectrum 
hazardous also in the infrared. Further, as the effective tem- 
perature decreases the peak of the spectra moves redward, 
until at T e ff ~ 3000 K it settles in J band and stops mov- 
ing further to the red (Allard & Hauschildt 1995). Below 
~ 4000 K, depending on the infrared band, Rtheo flattens 
out and then monotonically increases with decreasing effec- 
tive temperature. This behaviour in the infrared resemble 
that shown by Rtheo also in the optical bands (Figure [4}. 

It seems therefore obvious that below 4000 K, depending 
on the metallicity, particular care must be used in determin- 
ing the effective temperature by means of the flux ratio. On 
the other hand, since in cool stars Rtheo behaves qualita- 
tively the same in both infrared and optical colours, once a 
technique for determining T c s is found, that can be readily 
applied to any photometric band. 

As we discuss in more detail in Appendix A, when 
Rtheo monotonically increases with decreasing temperature 
it is still possible to converge in T e a if we compare the ob- 
served flux product TboI (Earth) x T\ (Earth) to its theoret- 
ical counterpart (|) <jT c ?^T\ (model). The apparent draw- 
back of this method is that it introduces a dependence on 
the angular diameter, whereas such a dependence was can- 
celed out when doing the flux ratio. However at each n th 
iteration the angular diameter can be estimated from the 



(n - 1) iteration 




Figure 4. Ratio between the bolometric and the monochromatic flux (Rtheo) in various filter bandpasses £ (lines of different colours) 
as function of T e g for the Phoenix models in three different metallicities +0.0 (solid), —1.0 (dashed) and —2.0 (dotted). Broad-band 
(heterochromatic) fluxes have been reduced to monochromatic fluxes as described in the Appendix B of Casagrande et al. (2006). 



T Bo i (Earth) (ri 



(2) 



so that it is still possible to converge in effective tem- 
perature 



^(Earth^^Tjtt^^ 
T Bo i (Earth) , _ 1} T\ (model) 



l(n-l) 



(3) 



Since this formalism is valid when Rtheo monotonically 
increases with decreasing temperature, the advantage of this 
approach is that now it is possible to use also the optical 
colours to converge in effective temperature below ~ 4000 K. 
Notice that to bootstrap either the IRFM or the MOITE 
one needs to interpolate over a grid of synthetic spectra 
according to the details given in Appendix A. To do so, 
the metallicity of a star must be known: when this is not 
possible, [M/H] has been obtained with the procedure we 
describe in the next Section. 



4.3 Estimating the metallicities of M dwarfs with 
the MOITE 

Figure [4] shows that going to cooler T e ff, both optical and 
infrared colours start to show a strong dependence on the 
metallicity. We exploit this particularity to implement a 
novel technique to estimate the metallicities of the M dwarfs. 

The method works as follows: for a given star of un- 
known metallicity, we apply the MOITE to recover the ef- 
fective temperature, assigning each time a trial [M/H] to the 
star, from —2.1 dex to 0.4 dex, in steps of 0.1 dex. The trial 
metallicity assigned to the star is used for interpolating over 
the grid of model atmospheres. The chosen metallicity range 
well brackets our (rather local) sample of M dwarfs. For a 



given star we obtain 6 which estimates T e ff from each of 
the colour bands individually V(RI)cJH Ks , for each of the 
26 different metallicity choices (from —2.1 to 0.4 dex). Since 
each band has a different sensitivity to the metallicity, the 
scatter among the 6 is at a minimum when the correct 
metallicity is chosen, as we prove in the next Section for a 
set of synthetic colours. 

For real data, there is an additional complication. Em- 
pirically, for the 118 M dwarfs with known [M/H], we find 
that the temperature estimates in each band Tj are on av- 
erage offset by a few 10s of K from the average T e g with a 
dependence on both [M/H] and T g. This offset might be 
ascribed to zero-point errors in the absolute calibration (see 
also fig. 8 and 12 in Casagrande et al. 2006) and/or to sys- 
tematics in the spectral library, or both. The computation of 
TboI (Earth) and T\ (Earth) from the observed multi-band 
photometry depend on the adopted zero-points and absolute 
calibration (Casagrande et al. 2006). Although our adopted 
Vega absolute calibration has been thoroughly tested in both 
optical and infrared bands via ground (Tokunaga & Vacca 
2005) and spaced based (Bohlin & Gilliland 2004; Price et 
al. 2004; Bohlin 2007) measurements, uncertainties at the 
level of a few percent are present and are almost certainly 
responsible for the systematic offsets of order tens of Kelvin 
in between different bands. Stars with different metal- 
licity and effective temperature emit differently in a given £ 
band; since the adopted absolute calibration and zero-points 
change the contribution of each £ band into the final result, 
this explains why the offsets are function of [M/H] and T c ff . 
We use the 118 M dwarfs for which we know their metallic- 
ities to correct these offsets in obtained from each band. 
For the real stars, this reduces the scatter in temperatures 
for each of the trial metallicities, and considerably assists 
in the recovery of the correct metallicity. Notice that this 
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Figure 5. Upper panels: accuracy of the MOITE in recovering the correct [M/H] and T c ff when realistic observational uncertainties 
(Section 14.41 1 are included in the Monte Carlo simulation. Different line-styles correspond to the recovered parameters for T e ff ^ 4300 K 
(continuous line), 3900 K (dotted) and 3500 K (dashed). Lower panels: contour plots showing the accuracy of the method as function of 
T e ff . When going to cooler effective temperature the fluxes in different bands become extremely sensitive to the metal content and the 
physical parameters of the underlying synthetic spectra are always recovered with excellent accuracy, as discussed in the text. 



correction to properly estimate [M/H] for our stars is cali- 
brated on the Bonfils et al. (2005) metallicity scale, but for 
a given metallicity, T e g is obtained with the MOITE alone. 



4.4 MOITE : accuracy of the technique 

The first test for the MOITE is to ensure that the proposed 
technique works, and if so, to which accuracy. The best way 
to address the level of internal accuracy of the method is by 
using synthetic colours to check whether the correct physical 
parameters (i.e. T e g and [M/H]) of the underlying synthetic 
spectra are recovered. 

We use the Phoenix model atmospheres to generate a 
set of synthetic BV (RI)c J H K s magnitudes for stars in the 
temperature range 2100 ^ T e g ^ 4500 K and metallicity 
range —2.0 ^ [M/H] 0.0 dex. Notice that now the adopted 
absolute calibration and zero-points are not responsible for 
any offset among T| in different bands, since the same abso- 
lute calibration and zero-points are used in generating syn- 
thetic magnitudes and in the MOITE. 

We begin by applying the MOITE, regarding [M/H] as 
a fixed, known parameter. To get the iteration started, ini- 
tial temperature estimates were made using Bessell's (1991) 
Teg : (R — I)c calibration. We find that we can recover 
the correct effective temperatures of the model spectra with 



an accuracy of 1 ± 3 K and the bolometric luminosities (i.e. 
aT e fi 4 ) within 0.1 percent. We then tested what happens if 
rather than using the Bessell (1991) T e g : (R — I)c calibra- 
tion for the first estimate of the effective temperature we 
start from either T cS = 5000 K or from 2000 K. The method 
still recovers the correct temperatures and bolometric lumi- 
nosities with the same accuracy as before (but more itera- 
tions are needed). The MOITE is thus pretty insensitive to 
poor first guesses of the effective temperature and always 
correctly converges. 

We have then tested the MOITE at recovering metal- 
licity (i.e. introducing [M/H] as a free parameter) as well as 
effective temperature and luminosity. We find we can recover 
the metallicities of the underlying synthetic spectra with an 
accuracy of 0.006 ±0.04 dex, with very few cases when they 
deviate by 0.1 dex. As a consequence, T e ff and the bolomet- 
ric luminosities are still recovered with very good accuracy. 
These tests establish that the technique has a high internal 
accuracy in the ideal case of no observational errors. 

For a more realistic approach, we mimic real data by 
running Monte Carlo simulations with realistic observa- 
tional uncertainties. For a set of synthetic BV (RI)c J H K s 
colours, we have assigned each time random errors with a 
normal distribution centered on the synthetic values and a 
standard deviation equal to the typical optical (Section l2.ll) 
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Table 1. Observable and physical quantities for our sample stars. 



name T eff ± AT cff 9 ± A(9 m Bol V B-V U-B V - R c V - I c J H K s [M/H] 

(K) (mas) 

HIP112 3923 ± 142 0.209 ±0.016 9.668 10.748 1.410 1.320 0.879 1.722 8.017 7.408 7.217 -0.08 

HIP897 3786 ± 126 0.228 ±0.016 9.628 10.822 1.463 1.249 0.914 1.826 7.976 7.314 7.119 -0.18 

HIP1734 3397 ± 154 0.306 ± 0.029 9.459 11.133 1.508 1.162 1.009 2.209 7.674 7.052 6.785 0.07 

HIP1842 3327 ± 110 0.241 ±0.017 10.073 11.886 1.523 1.166 1.042 2.326 8.259 7.640 7.375 0.04 

Effective temperatures (T e fj) and angular diameters (8) are computed as described in the text. For the stars with HIP number, the 
metallicities are from the Bonfils et al. (2005) calibration, while the remainder are obtained with the MOITE. We only give metallicitics 
for stars with T e g above 3080 K for the reasons explained in Section [5] Apparent bolometric magnitudes (ms ;) are obtained according 
to Section [6] where the absolute bolometric magnitude of the Sun MboI q = 4.74. Optical colours are in the Johnson-Cousins system, 
whereas infrared are from 2MASS. The full table is available in electronic format. 



and infrared (Section I2.2[l photometric errors. The results 
are shown in Figure[5] and demonstrate that the method be- 
comes more accurate when going to lower T e g, since at cool 
temperatures the flux ratios show a pronounced dependence 
on the metallicity as expected from Figure [4] Now that real- 
istic observational errors are taken into account we recover 
the metallicities of the synthetic spectra within 0.1 — 0.3 dex 
and the effective temperatures within 50 — 100 K, also de- 
pending on the spectral type. It is important to remember 
that real data might include systematic uncertainties (espe- 
cially in the absolute calibration) which are difficult to assess 
and therefore not included when running the Monte Carlo 
simulations. The actual accuracy might be somewhat worse 
than in Figure [5] Also, it is fair to remember that the com- 
parison with the observed colours (Section 13 - 1 ft shows that 
at coolest temperatures there is still room for improvements 
in theoretical modeling. For very cool stars, our results are 
therefore subject to possible refinements. However, even in 
the range 3000 < T cS < 4000 K the method recovers [M/H] 
and T e ff of the underlying synthetic spectra with an accu- 
racy of ~ 0.2 dex and ~ 50K at a confidence level of 85 
percent (i.e. within 1.5 a). 

We conclude that the MOITE shows a high degree of 
internal accuracy over the range of metallicities and effective 
temperatures covered in the present study. However, for a 
better estimate of the reliability of such a technique, we 
need to compare our results with direct measurements. This 
will be done in Section and [7J where we will compare the 
metallicities and the effective temperatures returned by the 
MOITE with those recently measured in literature by means 
of other techniques. 

4.5 The final error budget 

We have shown our technique to be very promising in ob- 
taining both effective temperatures and metallicities of cool 
stars. Based on the comparison with other measurements, 
in the next Section we will estimate the uncertainty of our 
metallicities to be on average 0.2 — 0.3 dex especially when 
photometry is available in many bands. 

For evaluating the final errors in effective temperatures, 
bolometric luminosities and angular diameters, we use the 
same prescriptions as in Casagrande et al. (2006). Briefly, for 
each star we run 200 Monte Carlo simulations assigning each 
time random errors to the photometry (according to the un- 
certainties given in Section [2.11 and | 2,2[) and to [M/H]. We 



have also accounted for a change of ±0.5 dex in the value of 
log(g) used in our grid of model atmospheres and for the case 
the errors in the absolute calibration correlate to give sys- 
tematically higher or lower fluxes. When using the MOITE 
to estimate the metallicities, we search for the solution that 
minimize the scatter among T%. Although the scatter is min- 
imized, it is always finite: this reflects photometric errors as 
well as possible systematics which are not fully corrected, 
especially at the lowest temperatures. We adopt a very con- 
servative approach and also include the scatter in different 
bands in the final T c g uncertainty. 

The effective temperatures, msoi and the metallicities 
for our entire sample are given in Table [1] On average, our 
effective temperatures are accurate to about 100 K which 
corresponds to 2 — 3 percent in the studied temperature 
range. The uncertainty in bolometric luminosities is usually 
between 3 and 4 percent and that in angular diameters be- 
tween 4 and 6 percent. Also, the fact that the MOITE recov- 
ers the metallicities consistently with other determinations, 
although is not a proof of the correctness of T e g is reassur- 
ing and indirectly confirms the accuracy of the temperatures 
(Figure [5]). 



5 THE M DWARF METALLICITY SCALE 

In the previous Section we have introduced how the MOITE 
works and we have evaluated its accuracy by means of Monte 
Carlo simulations. Here, we compare the metallicities ob- 
tained with our technique to those recently measured by 
Bonfils et al. (2005), Woolf & Wallerstein (2005, 2006) and 
Bean et al. (2006a, 2006b). For dwarfs with T cff below 3000 K 
we use Hipparcos common proper-motion companions to 
compare the metallicity we derive for the cool secondaries 
with that more easily measured for the hotter primaries in 
the aforementioned studies. We will compare our effective 
temperatures with other determinations available in the lit- 
erature in Section [7J 



5.1 Accuracy of the metallicities above 3000 K 

Figure [6] shows our metallicity estimates for the 118 M 
dwarfs with known metallicities, and shows a la scatter of 
0.2 dex, i.e. within the accuracy of the Bonfils et al. (2005) 
calibration. This comparison is encouraging, but still in- 
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the one-to— one relation. Lower panels: the metallicity difference 
([M/H] Bonfila — [M/H] MOITE ) as function of other parameters. 
Long-dashed lines are the la scatter. 



eludes stars which were used to construct the calibration 
(Section l4.3[) . so is not a fully external check on the method. 

Recently, reliable metallicities for M dwarfs have been 
measured by Woolf & Wallerstein (2005, 2006) from very 
high-resolution spectra. Spectral synthesis technique has 
been also successfully applied by Bean et al. (2006a, 2006b) . 
We have extensively searched for BV (RI)c J H K s photom- 
etry of M dwarfs analyzed in the aforementioned studies and 
found accurate measurements for those reported in Tabled 
for T e ff approximately above 3000 K the mean difference in 
metallicity is just — 0.03±0.06 dex (a — 0.17 dex). In partic- 
ular, the comparison with the direct spectroscopic measure- 
ments of Woolf & Wallerstein (2005, 2006) agrees always 
within 0.13 dex (filled diamonds in Figure [7J). The agree- 
ment with the [M/H] measurements of Bean et al. (2006a, 
2006b) is somewhat poorer (asterisks in Figure [7J, but there 
are large discrepancies between their and our T e fr scale, as 
we discuss later in this Section. Such large differences in the 
adopted T c s obviously reflect in the abundances measured 
by Bean et al (2006a, 2006b). 

Notice that our technique relies on the metallicities de- 
termined from Bonfils et al. (2005), whose calibration also 
includes measurements from Woolf & Wallerstein (2005). 
To further and independently test our results, we have ap- 



plied the MOITE to all the stars in Woolf & Wallerstein 
(2006). For these stars, only VJHKs photometry is avail- 
able. We have used V magnitudes as given in table 1 of Woolf 
& Wallerstein (2006) and JHK S from 2MASS. Again, we 
have used only stars with total photometric errors in the in- 
frared smaller than 0.10 mag. Since we are now running the 
MOITE using fewer colours, we might expect our results to 
be slightly less accurate. The comparison with the Woolf & 
Wallerstein (2006) measurements in Figure [7] is reassuring, 
with a mean difference of 0.20 ±0.05 dex (a = 0.27 dex) and 
it validates our technique over a larger range of metallicities. 

Differently from the study of FGK dwarfs, M dwarfs 
exhibit complex spectra which render far less trivial the de- 
termination of both effective temperatures and metallicities 
from spectroscopic analysis only. The purpose of this Sec- 
tion is to evaluate the accuracy of our metallicities, whereas 
we compare our temperature scale with other existing ones 
in Section [7] It is however worth discussing here the differ- 
ences with the effective temperatures adopted by Woolf & 
Wallerstein (2005, 2006) and Bean et al. (2006a, 2006b) to 
derive their metallicities. Our effective temperatures are sys- 
tematically cooler than those of Woolf & Wallerstein (2006) , 
especially for earlier M spectral types (3700 K and above), 
where the difference in metallicities is also higher. Woolf & 
Wallerstein (2005, 2006) estimate T e s from theoretical color- 
temperature relations in V — H and V — Ks obtained using 
older model atmospheres (Hauschildt et al. 1999). Accord- 
ing to Woolf & Wallerstein (2006), systematics as high as 
100 — 200 K in their T c a could not be excluded. Our effec- 
tive temperatures are in agreement with the latest Phoenix 
and above 3500 K Castelli & Kurucz (2003) models, as it 
can be seen from Figure [9] The disagreement in B — V for 
the Castelli & Kurucz (2003) models is likely due to pitfall 
in accounting the contribution of all molecular features in 
these spectral regions. The effective temperatures obtained 
by the spectral fitting technique of Bean et al. (2006a, 2006b) 
rely particularly on the TiO bandhead. Bean et al. (2006a) 
also notice that while at M0. 5 spectral types their T e ff agree 
with the scale of Reid & Hawley (2005), their effective tem- 
peratures increase linearly with spectral type and for their 
latest M dwarf (GJ 105 B) they are about 300 K hotter than 
Reid & Hawley (2005) and 400 K hotter then our scale. This 
hotter temperature scale is possibly related to drawback in 
determining effective temperatures via spectral synthesis of 
certain bandhead (e.g. Jones et al. 2005). We will further 
discuss and test our temperature scale with other recent 
determinations, in particular with interferometric angular 
diameters, in Section[7]and prove our effective temperatures 
to be reliable. 

Concluding, the metallicities estimate with the MOITE 
above ~ 3000 K agree within 0.2 — 0.3 dex to those measured 
by other recent studies, especially when all BV(RI)cJHKs 
photometry is available. In particular, the comparison with 
direct spectroscopic measurements of Woolf & Wallerstein 
(2005) suggest an even better agreement, at the level of 
0.1 dex. 



5.2 Accuracy of the metallicities below 3000 K 

As we have already mentioned, the calibration performed in 
Section f4. 31 is function of T a g and [M/H] and in our case is 
obtained in the metallicity and temperature range shown in 
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Table 2. Magnitudes and colours for various M dwarfs and the T c g and [M/H] recovered with the MOITE compared with those obtained 
by other techniques. In the infrared we have used JHK$ magnitudes from 2MASS. 



Name NLTT V B-V V - R c V - I c Ref. T cff MOITE (K) [M/H] MOITE T e g (K) [M/H] mcas . Ref. 



GJ 191 


14668 


8.841 


1.570 


0.956 


1.951 


1 


3661 


-1.00 


3570 


-0.99 


a 


GJ 701 


45883 


9.362 


1.515 


0.976 


2.060 


2 


3557 


-0.20 


3630 


-0.20 


a 


GJ 828.2 


51282 


11.090 


1.534 


0.964 


1.967 


2 


3646 


-0.24 


3680 


-0.37 


b 


GJ 876 


55130 


10.179 


1.571 


1.182 


2.733 


2 


3076 


-0.09 


3478 


-0.12 


c 


GJ 581 


39886 


10.568 


1.602 


1.109 


2.501 


2 


3211 


0.00 


3480 


-0.33 


c 


GJ 297.2 B 


19072 


11.80 


1.49 


1.03 


2.29 


4 


3370 


-0.20 


3659 


-0.12 


cl 


GJ 105 B 


8455 


11.66 


1.50 


1.22 


2.78 


4 


3048 


-0.30 


3444 


-0.09 


cl 


GJ 412 B 


26247 


14.44 


2.08 


1.66 


3.77 


3 


2700 


-0.31 




-0.43 


a 


GJ 618 B 


42494 


14.15 


1.79 


1.412 


3.233 


3 


2614 


-0.26 




-0.12 


e 


GJ 752 B 


47621 


17.20 




2.10 


4.36 


3 


2250 


-0.04 




-0.04 


c 



Source of the optical photometry : (1) Kilkenny et al. (1998); (2) Koen et al. (2002); (3) Bessell (1990a); (4) Laing (1989). Source of 
metallicity : (a) Woolf & Wallerstein (2005); (b) Woolf & Wallerstein (2006); (c) Bean et al. (2006b); (d) Bean et al. (2006a); (e) Bonfils 
et al. (2005). For the Bonfils et al. (2005) calibration we have used the V magnitudes given in the Table and the Kg magnitudes from 
2MASS. 



Figure [6] The comparison in Figure [7] suggests that the tech- 
nique can be safely applied to metallicities as low as about 
—2.0 dex, at least above 3500 K and to metallicities typical 
of the solar neighbourhood down to about 3000 K. Here, we 
would like to test to what extent we can use the MOITE 
to measure metallicities in stars with effective temperatures 
below this latter limit. The extrapolation to lower temper- 
atures is particularly interesting, since on the basis of our 
Monte Carlo simulations the MOITE is expected to recover 



the metallicity to high accuracy (always within 0.1 — 0.2 dex) 
below 3000 K (Figure [5}. Also, for the coolest effective tem- 
peratures, as we discussed in Section T3. II theoretical models 
do show discrepancies with respect to the empirical colour- 
colour diagrams. Those could reduce the formal high accu- 
racy expected for T e a < 3000 K, so that below this temper- 
ature it is interesting to compare our with respect to other 
more "direct" metallicity estimates. 

We are not aware of any [M/H] measurement for M 
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Figure 8. Upper and central panel: bolometric corrections in V and R bands as function of various color indices for our M dwarfs. 
Lower panel: relation between the magnitudes observed in different bands and the bolometric magnitude mg [. Linear fits in the 
form mg i = A + B are given along with the resulting standard deviation a. 



dwarfs with T cff < 3000 K, but from the list of Gould & 
Chaname (2004), we have searched for multi-band optical 
and infrared photometry of late M dwarfs which are com- 
mon proper-motion companions to Hipparcos stars. We have 
then compared our metallicities estimated with the MOITE 
to those more easily measured for the primary Hipparcos 
stars. We have found three late M dwarfs (GJ 412 B, GJ 
618 B, GJ 752 B) whose common proper-motion compan- 
ions are early type M dwarfs with metallicities from Woolf 
& Wallerstein (2005) or from the calibration of Bonfils et 
al. (2005) . We caution that some of the stars in Table [2] are 
classified as "variable" in SIMBAD and GJ 618 B has very 
large photometric errors in 2MASS (much larger than our 
usual fiducial level "h_" and "k_msigcom" < 0.10 mag). 

There are only three stars and it is hard to draw defini- 
tive conclusions, but it is reassuring to see such a nice agree- 
ment, with the metallicity of the secondaries in agreement 
with that of the primaries within 0.14 dex even in the worst 
case. For two stars [M/H] of the primary is obtained using 
the Bonfils et al. (2005) formula, which we have used to cal- 



ibrate our technique. Nonetheless, now we are working with 
effective temperatures much cooler than those of Figure [6] 
and the metallicities are still properly recovered. Of course 
more stars are needed before confidently extend our tech- 
nique to such cool stars, but the method looks promising. 
To be on the safe side, for our entire sample of M dwarfs (Ta- 
ble PJ) we give [M/H] only for stars with T cS above 3080 K, 
i.e. only where the MOITE is safely calibrated as explained 
in Section \4. 31 



6 EMPIRICAL BOLOMETRIC CORRECTIONS 

We adopt the same definition of Casagrande et al. (2006) to 
define the bolometric correction in a given £ band, where 



BC^ = rriBoi — m£ 



(4) 



and the zero-point of the msoi scale is fixed by choosing 
MboI,q = 4.74. Although it is possible to give analytic trans- 
formations between the flux and various colour indices, bolo- 
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metric corrections can be readily computed using equation 
([4]) in any band from the stars in Table [1] and are probably 
more useful. 

The bolometric correction in V and R bands as func- 
tion of various colour indices are shown in Figure [8] In the 
optical there is some dependence on the metallicity among 
the coolest stars, but when going to longer colour baselines, 
the data show very tight relations, especially in V — J. The 
flux emitted by cool stars peaks in the near-infrared. The 
bolometric corrections in Ic, J, H and Ks are almost con- 
stant as function of different color indices. For these bands 
it is therefore possible to pass directly from the observed 
to the bolometric magnitude m,Boi via linear fit as shown in 
Figure [8] and given here: 



m,Boi 



1.46 + 0.90 m Ia 
1.22 + 1.07 mj 
1.94 + 1.06 ms 
2.07 + 1.08 m Ks . 



(5) 



Of course, by using the data in Table [T] to fit the bolometric 
correction as function of a given colour index it is possi- 
ble to achieve higher accuracy still. Nonetheless, the linear 
fits in Figure [8] are useful to have a quick estimate of the 
bolometric magnitude given the apparent one. From these 
bolometric corrections and the colour-temperature relations 
given in the next Section it is also possible to obtain very 
accurate estimates of the angular diameter of M dwarfs, as 
we describe in Section [7. II 



7 THE M DWARFS TEMPERATURE SCALE 

We have fitted the observed colour — T c g relations by ana- 
lytical fits, using the following functional form: 



9 e « = Q.o + a-iX + a 2 X + a 3 X 



(6) 



where c g = 5040/T e gf, X represents the colour and en (i = 
0,1,2,3) are the coefficients of the fit. Depending on the 
band (Figure O, the data show an increasing scatter in 
the effective temperature estimates below ~ 2800 K, which 
is very likely due to metallicity dependencies. Because of 
this, we strongly caution against using V — Rc, V — Ic and 
(R — I)c indices below 2800 K. We have also verified that 
such scatter at low temperatures is not present in the other 
transformations of Table [3] We do not have enough stars be- 
low this temperature to address the question further and the 
metallicity estimates are still somewhat uncertain so that we 
do not include any metallicity dependence in our analytic 
fits. Furthermore our M dwarfs are drawn from the solar 
neighbourhood, and they have a limited range of metallic- 
ities. It will be interesting to try the method on halo M 
dwarfs, as these become available especially in large photo- 
metric/spectroscopic surveys currently underway (e.g. from 
RAVE, SEGUE, SkyMapper) and later with GAIA. We thus 
differ from the Casagrande et al. (2006) fitting formulae, 
which accounted for the larger metallicity coverage of that 
sample. Also, now we need to fit a third order polynomial 
to account for the inflection that occurs at lower T e g. For 
this reason, the fitting formulae of Table [3] are not an ex- 
act continuation of those in Casagrande et al. (2006), which 
are correct for T B g hotter than 4400 K. Therefore, the colour 
ranges for the fits in Table[3]do not overlap with Casagrande 



et al. (2006), but are given for slightly redder colours. If a 
link between the two scales is needed, however, we advise 
the users to a careful case by case study. 

With the exception of B — V, none of the colour- 
temperature transformations have strong dependence on 
[M/H] above 2800 K, and therefore our relations are not 
likely to be affected much as metallicities for M dwarfs im- 
prove. The temperature fit to the B — V colour has huge 
scatter, and gives only the crudest temperature estimate. If 
one really needs to estimate T c g from B — V the best choice 
is probably to use Figure 

We have searched for DENIS photometry so as to give 
in Table [3] the colour — T e ff relations also in this system, 
although for a smaller number of stars. According to the 
DENIS database, those magnitudes have larger photometric 
errors (on average between 0.05 to 0.09 magnitudes) than 
2MASS so this might explain why the relations in Table [3] 
are less accurate for the DENIS colours. Also, we have found 
relations involving the /denis very noisy and we do not give 
them. 

For most of the stars in Koen et al. (2002) SAAO JHK 
photometry is also available from Kilkenny et al. (2007) and 
in Table [3] colour-temperature relations are given in this 
system, too. Stars with SAAO infrared photometry are all 
hotter than 3000 K and because of the reduced temperature 
range, second order polynomial fits are accurate enough now. 
We caution against extrapolating these relations to cooler 
temperatures, in particular for indices Ic — Jsaao, Ic — 
-Hsaao and Ic — -Ksaao 

Many M dwarfs are intrinsically variable, owing to spots 
or other activity. We have cleaned the sample from variable 
stars as best we can, (and indeed our colour-temperature 
relations are very tight), however there might still be unrec- 
ognized long term variables present. The relations we give 
thus apply also to intrinsically low variability stars. 

For the sake of completeness, we have applied the 
MOITE to the M dwarfs labeled as variable in the original 
Koen et al. (2002) sample and which were excluded, accord- 
ing to the selection criteria of Section]!?] The stars were found 
to follow the same mean loeii of the colour-temperature rela- 
tions as the non- variables, but with a larger scatter, usually 
about twice the a(T e g) of TableO One should keep in mind 
that the stars in Koen et al. (2002) were already preselected 
in Hipparcos in order to exclude extremely active stars: how 
well our results would apply to these latter objects we leave 
to future studies. 

The latest Phoenix model atmospheres show good 
agreement with our colour-temperature relations in Figure 
[9] This partly reflects the fact that our scale has been ob- 
tained using these models in the MOITE, but we are using a 
great deal of observational information to recover the total 
bolometric luminosity, and the model atmospheres are used 
only to estimate the missing flux, which is at most of order 
twenty percent. We demonstrate further the reliability of our 
temperature scale in what follows: in Section 17.11 using re- 
cent determinations of angular diameters for M dwarfs and 
in Section \7. 2 1 comparing our results with those obtained by 
various recent temperature studies. 
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Table 3. Coefficients and range of applicability of our colour-temperature relations (Eq. [B). The photometric systems are V(RI)c for 
Johnson-Cousins, JHKg for 2MASS, (JA)denis for DENIS and (JHK)saao for SAAO. For some indices, we caution the users from 
extrapolating these relations to redder colours than those given, as explained in Section [7] 
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Notes. N is the number of stars employed for the fit after the 3<r clipping and cr(T c ff ) is the final standard deviation (in Kelvin) of the 
proposed calibrations. 



7.1 MOITE versus interferometric angular 
diameter measurements 

Although in the past, much work has been done in determin- 
ing the effective temperature scale of the M dwarfs, a firmly 
established scale has not been achieved. Until recently, in 
fact, the only two M dwarfs with measured linear diameters 
were the eclipsing binaries YY Gem (Kron 1952, Habets & 
Heintze 1981) and CM Dra (Lacy 1977, Metcalfe et al. 1996), 
but the limiting factor in accurately determining their effec- 
tive temperatures were the parallaxes. Long-baseline inter- 
ferometry has recently provided angular diameters measure- 
ments for an handful of nearby M dwarfs which can be used 
to test the accuracy of our effective temperature and bolo- 
metric luminosity scale (Segransan et al. 2003; Berger et al. 
2006). Unlike G and K dwarfs, for which all the interfero- 
metric targets have saturated 2MASS photometry, half of 
the M dwarfs with measured angular diameters have good 
2MASS colours. 

The stellar angular diameters obtained with the 
MOITE are computed from the basic definition 



TboI (Earth) = 




so that in principle a conspiracy of wrong effective temper- 
atures and bolometric luminosities could still return correct 
angular diameters. However, the bolometric luminosities of 
our targets are observed via multi-band photometry (only 
subject to minor corrections to estimate the missing flux, 
see also Appendix A), so that TboI (Earth) is fixed by the 



observations and therefore comparison of our angular diam- 
eters with those measured by interferometers automatically 
tests our temperature scale. 

We caution that even interferometric angular diame- 
ter measurements depend mildly on modelling assumptions, 
in particular the limb-darkening corrections to convert the 
measured uniform-disk angular diameters into the physical 
limb-darkened disks (#ld) and to which we compare our 8 
of equation ©. The limb-darkening coefficients used for M 
dwarfs (Claret 2000) are computed using solar abundance 
atmospheric models, whereas the interferometric targets of 
Table [4] span a larger metallicity range. Another source of 
uncertainty is due to the fact that limb-darkening coeffi- 
cients are calculated using ID atmospheric models, whereas 
3D models predict a less significant center-to-limb variation. 
Such a difference might be up to a few percent in #ld for 
hotter F and G stars, but is expected to be much smaller in 
the case of M dwarfs (Allende Prieto et al. 2002; Bigot et al. 
2006). Since all interferometric measurements reported here 
have been performed in the infrared, where limb-darkening 
effects are minimized, we expect these uncertainties to be 
within the observational errors. 

When running the MOITE for M dwarfs in Table g] 
with good 2MASS photometry we have used the metallici- 
ties from Woolf & Wallerstein (2005) or applied the Bonfils 
et al. (2005) calibration, except for GJ 699 which is out- 
side of the Bonfils' et al. (2005) range of applicability and 
for which a solar metallicity is thought to be appropriate 
(Leggett et al. 2000; Dawson & De Robertis 2004). The er- 
rors have been computed as described in Section 14.51 For 
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Figure 9. Colour — T c ff plots in different bands for our M dwarfs. Overplotted are the prediction from the Phoenix models (yellow lines) 
for two different metallicities which roughly bracket our sample of stars. Also shown for comparison the prediction from the Castclli & 
Kurucz (2003) models for solar metallicity (blue line). Red squares in the T e g vs. V — Ic plot are from the temperature scale of Reid & 
Hawley (2005). 



Table 4. Comparison between the MOITE effective temperatures and angular diameters (columns 8 and 9) and the interferometric 
measured ones (columns 10 and 11). 



Name 


V 


U -B 


B-V 


V-Rc 


V-I C 


Ref. 


T cff (K) 


0MOITE 


T cS (K) 


6>LD 


Ref. 


GJ 


191 


8.841 


1.186 


1.570 


0.956 


1.951 


1 


3661 ± 77 


0.637 ±0.028 


3570 ± 156 
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s 


GJ 
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GJ 


411 f 
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2.15 
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3467 ± 104 
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3570 ± 42 


1.436 ± 0.030 
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GJ 
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9.05 




1.52 


0.98 


2.04 
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3594 ± 101 


0.636 ± 0.037 


3243 ± 160 


0.753 ± 0.052 
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GJ 


526 f 


8.464 
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2.070 


5 


3533 ± 106 


0.884 ± 0.053 


3636 ± 163 


0.845 ± 0.057 
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GJ 
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3145 ± 69 
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4 


3544 ± 153 
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3577 ± 107 
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3626 ± 56 
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s 



Source of the optical photometry: (1) Kilkenny et al. (1998); (2) Koen et al. (2002); (3) Celis (1986); (4) Bessell (1990a); (5) The, 
Steenman & Alcaino (1984); (6) Landolt (1983). In the infrared we have used 2MASS JHKg photometry (not reported here). A f 
indicate a poor 2MASS photometry, so that the angular diameter has been obtained from the calibration of Section [JJ and [6] and not 
running the MOITE directly. Source of interferometric measurements: (s) Segransan et al. (2003); (b) Berger et al. (2006). 



the other stars (i.e. those with inaccurate 2MASS photom- 
etry) we apply the bolometric luminosity and effective tem- 
perature calibrations of Section [6] and [7J We estimate T e s 
using the V — Ic index, which has little intrinsic scatter 
above 2800 K (see Tableland Figure 0. We then compute 
the bolometric correction in V band by linearly fitting the 
BCv vs. V — Ic relation (Figure[8]) in the colour range [1.95, 
2.25]. As for the temperature calibration, this relation has 



little intrinsic scatter and the linear fit in the given range 
is accurate to 0.015 mag. Once the bolometric correction 
and the effective temperature are known, the angular diam- 
eter can be readily computed (eq. 14 in Casagrande et al. 
2006). To these stars, we assign a 3 percent error in T e ff and 
6 percent error in angular diameter, consistently with the 
errors obtained for the other stars in the sample. Although 
for these stars we are not applying the MOITE directly, we 
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are using the calibrations obtained with the MOITE itself 
so that they are fully representative of our temperature and 
luminosity scale. 

The overall comparison between our and the interfero- 
metric angular diameters shown in Figure [10] is very good. 
Segransan et al. (2003) and Berger et al. (2006) also com- 
pute effective temperatures which are obtained combining 
the measured #ld with bolometric correction polynomial fit 
or comparing observed and model-predicted fluxes based on 
the observed angular diameters. This allows a direct com- 
parison with their effective temperatures in Figure [TT1 (filled 
diamonds) . Only two stars deviate by more than la from the 
one-to-one relation in Figure [TUI and [TT1 namely GJ 514 and 
GJ 880. For these two stars, however, the two methods of 
determining T c ft from 6>ld (i.e. using bolometric correction 
or comparing observed and model-predicted fluxes) return 
results discordant by 100 K. Choosing T c g computed using 
bolometric correction (table 4 in Berger et al. 2006) would 
reduce by 100 K the discrepancy for these two stars in Fig- 
ure [11] We do not know the reason for such disagreement, 
however, our bolometric corrections are carefully determined 
from multi-band data, whereas the main point of interfer- 
ometric studies is to precisely measure angular diameters 
(with which we are in good agreement), rather than deter- 
mining accurate empirical bolometric corrections. The im- 
portant point from the comparison in this Section is that 
overall our angular diameters and effective temperatures 
agree with the trend defined by interferometric studies. 

In particular, Barnard's star (GJ 699) is one of the 
benchmarks in setting the cool star temperature scale. Our 
angular diameter is in excellent agreement with the interfer- 
ometric one and our effective temperature Toft = 3145 ± 69 
also closely matches the value of 3134 ± 102, obtained by the 
careful analysis of Dawson & De Robertis (2004) . 

7.2 Comparison with other temperature scales 

The determination of effective temperatures by means of 
different techniques below 4000 K becomes increasingly dif- 
ficult as a result of the increasing complexity of the stellar 
spectra and previous studies have usually computed T e s for 
a limited number of stars. Interferometry has recently pro- 
vided a breakthrough to anchor the temperature scale down 
to 3000 K, as discussed in Section r7.ll Moreover, stars cooler 
than 3000 K have angular diameters too small to be resolved 
by currently available interferometers. Here we compare our 
effective temperatures to those obtained by various recent 
studies. 

Among the possible different techniques to estimate 
T c ft , one is to fit observed molecular features with model 
predictions (e.g. Kirkpatrick et al. 1993). However, these fea- 
tures might depend not only on the effective temperature, 
but on other effects of line formation and the reliability of 
the models themselves. It is a long-standing result that espe- 
cially below 3500 — 3000 K such a technique returns effective 
temperatures that are higher by several hundred kelvin with 
respect to other more empirically motivated methods (e.g. 
Pavlenko & Jones 2002; Jones et al. 2005). Very recently Viti 
et al. (2008) have proposed a new and promising technique 
based on very high resolution mid-infrared observations of 
pure rotational water vapour transitions in M dwarfs. An- 
other approach is to compare observed and synthetic spectra 
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Figure 10. Comparison between our and the interferometric 
$LD angular diameters. Circles refer to angular diameters ob- 
tained applying the MOITE directly. Squares are for those stars 
for which the angular diameters have been computed from the 
colour— temperature and colour-luminosity calibrations of Section 
[7] and \E\ Dotted lines are intended to guide the eye. 



and estimate T e g from the model that better matches the 
observation in the infrared (e.g. Leggett et al. 2000, 2001), in 
the optical (e.g. Dawson & De Robertis 2000) or throughout 
most of the spectrum (e.g. Pavlenko et al. 2006). Burgasser 
& Kirkpatrick (2006) have shown that the parameters de- 
rived using optical or near-infrared fits for a given object 
exhibit clear systematic differences up to 100 — 200 K in ef- 
fective temperature and 0.5 — 1.0 dex in metallicity. 

A more consistent way to determine T e ft is to analyze 
the entire spectral region contributing to the bolometric flux, 
although in the past this approach has been done mostly 
with black-body calibrations rather than with model pre- 
dictions (e.g. Veeder 1974; Reid & Gilmore 1984). A more 
rigorous attempt to recover the bolometric flux was used by 
Tsuji et al. (1996). Finally, Reid & Hawley (2005) have col- 
lected spectroscopic and photometric T e g estimates of a few 
well studied nearby M dwarfs covering the spectral type M0 
to M9 (their table 4.1). 

We have searched in the literature for other T c g deter- 
mination of our stars and did not find many, particularly 
below 3000 K. To increase the number of stars available for 
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Figure 11. Comparison between the MOITE effective temper- 
atures and those determined in other recent studies. A typical 
error bar of ±100 K is shown in the lower right corner. Dotted 
line with slope one is intended to guide the eye. 



the comparison, we have applied our technique to a few 
more very red dwarfs (LHS36, LHS68, LHS292, LHS429, 
LHS2065, LHS2924, LHS3003) which are commonly studied 
in the literature. For these stars optical colours are available 
from Bessell (1990a, 1991) and infrared from 2MASS. We 
caution that they are all classified as flare stars in SIMBAD 
and for this reason they were not included in the original 
sample of Section [2] The effective temperatures for these 
and few other cool stars from Table [T] [5] and 3] are discussed 
below. 

LHS2 - (GJ 1002). For this star we obtain T cff = 2750 K, 
slightly cooler than the temperature of 2900 K obtained by 
Leggett et al. (2000). 

LHS36 - (GJ 406). We obtain an effective temperature 
of 2500 K which is cooler by about 100 K than the value of 
2600 K obtained by Leggett et al. (2000) and reported also 
in Reid & Hawley (2005). For the same star Pavlenko et al. 
(2006) obtained an effective temperature of 2800 K after a 
critical examination of the most recent model atmosphere 
fit to this object. The same temperature was also obtained 
by Tsuji et al. (1996). Golimowski et al. (2004) also found a 
hotter temperature than we do (2900 K). 

LHS37 - (GJ 411). For this star we obtain T cB = 3467 K 
which is in agreement with 3510 K reported in Tsuji et al. 
(1996), 3500 K in Leggett et al. (1996) and 3400 K in Reid 
& Hawley (2005). 

LHS39 - (GJ 412 B). We obtain an effective temperature 
of 2700 K, which again is cooler by 100 K with respect to the 
value obtained by Leggett et al. (2000) (2800 K). 

LHS57 - (GJ699). Our value T cff = 3145 K is halfway 
between 3100 K in Leggett et al. (2000), Reid & Hawley 
(2005) and 3210 K in Tsuji et al. (1996). 

LHS65 - (GJ 821). For this Hipparcos star we obtain 
?cff = 3567 K, in close agreement with 3600 K in Leggett et 
al. (2000). 



LHS68 - (GJ 866). This dwarf, together with GJ 406 is 
one of the reddest standards of the BV(RI)c system. Un- 
fortunately it is member of a triple system, which might de- 
crease the accuracy of the photometry (Delfosse et al. 1999). 
We obtain T c ff = 2650 K, considerably cooler than the value 
of 3000 K obtained by Dawson & De Robertis (2000) but 
in better agreement with T e ff = 2800 K in Reid & Hawley 
(2005). 

LHS292 - (GJ 3622). Our technique returns T cff = 
2450 K, cooler then both Leggett et al. (2000) (2600 K) and 
Golimowski et al. (2004) (2725 K). 

LHS429 - (GJ 644 C). The MOITE returns T cS = 2400 K 
for this late M dwarf, which is now 100 K hotter than in 
Leggett et al. (2000) but cooler than the value of 2500 K in 
Reid & Hawley (2005) and 2640 K in Tsuji et al. (1996). 

LHS473 - (GJ 752 A). The effective temperature we ob- 
tain 3343 K is midway between 3250 K in Reid & Hawley 
(2005) and 3475 K in Tsuji et al. (1996). 

LHS474 - (GJ 752 B). We obtain T cB = 2250 K which 
compares nicely with Tsuji et al. (1996) (2250 K) but it is 
slightly cooler than 2400 K in Reid & Hawley (2005). 

LHS2065 - (GJ 3517). For this very red dwarf we obtain 
T c ff = 2050 K in rough agreement with Leggett et al. (2001) 
(2100 K) but much cooler than the temperature of 2400 K 
obtained by Golimowski et al. (2004). 

LHS2924 - (GJ 3849). According to Reid & Hawley 
(2005) this is the best studied M9 dwarf for which they re- 
port T c ff = 2300 K. We obtain an effective temperature of 
2200 K which is slightly hotter than the value of 2130 K in 
Tsuji et al. (1996). 

LHS3003 - (GJ 3877). We obtain T cS = 2350 K which 
is now hotter than Leggett et al. (2001) (2200 K), but still 
cooler than Golimowski et al. (2004) (2600 K). 

Our temperatures are shown in Figure [TT] against other 
studies discussed here or in Section f7.1l There is an overall 
good agreement with respect to the effective temperatures 
obtained by Leggett et al. (2000, 2001), especially above 
3000 K. Below this temperature typical differences of order 
100 K exist, but on average the scatter in the data suggest 
we are on the same scale. Similarly, we are in very good 
agreement with the temperatures reported in Reid & Haw- 
ley (2005) above 3000 K, whereas below this temperature 
we are systematically cooler by 100 K. We also agree within 
100 K with Tsuji et al. (1996) except for effective temper- 
atures around 2500 K where we have two stars with larger 
differences. Our effective temperatures are also 250 — 400 K 
cooler than those in Golimowski et al. (2004) which are esti- 
mated using the relationship between the bolometric lumi- 
nosity (their observable) and Toff predicted from evolution- 
ary models. 

Summarizing, further work is needed before reaching a 
consensus among different temperature scales, especially be- 
low 3000 K. However, our temperatures are supported from 
interferometric angular diameters between 3100 and 3600 K 
and the homogeneous and smooth colour-temperature rela- 
tions of Figure [9] lead us to believe our temperature scale is 
credible below 3000 K. Also, the data in Figure [TT1 suggest 
that despite case-by-case differences, on average we agree 
with the effective temperature reported in many recent stud- 
ies. 
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8 ON THE DISCONTINUOUS TRANSITION 
FROM K TO M DWARFS 

In our previous paper we have implemented the IRFM to de- 
rive effective temperatures and bolometric iuminosities of G 
and K dwarfs. Here we have extended our technique to much 
cooler effective temperatures. For M dwarfs with accurate 
Hipparcos parallaxes (i.e. those from Koen et al. 2002) we 
can compute absolute magnitudes and therefore plot them 
with our previous sample from Casagrande et al. (2006) 
and study the properties of the entire lower main sequence 
(Figure IT2"j) . finding a very interesting feature: whereas the 
transition from late K to early M type occurs smoothly in 
the widely used My — (B — V) plane, in infrared colours 
(which are better tracers of Mb i and T e g) a prominent 
discontinuity appears around late K to early M types. The 
discontinuity is clearly evident in the observational plane 
Mk s — [V — K) (and we have verified it is present also in 
the other infrared colours) and is also quite prominent in 
the theoretical MboI — Tcft plane, at around 4200 — 4300 K 
(Figure HU). 

The feature is close to the point at which our two cal- 
ibrations (IRFM for G and K dwarfs and MOITE for M 
dwarfs) meet. To confirm the discontinuity is not depen- 
dent on the calibrations, we have reprocessed all the stars 
in Casagrande et al. (2006) with the MOITE, and using 
the Phoenix model atmosphere, confirming that we obtain 
the same temperatures and luminosities with both methods 
above ~ 4000 K. 

For both the GK dwarfs in Casagrande et al. (2006) 
and the M dwarfs studied here we have used strict selection 
criteria to remove double and variable stars. For the sake 
of completeness, we note that when the M dwarfs labeled 
as variable in Koen et al. (2002) are plotted in both the 
observational and Mb i — Tcff planes of the HR diagram, 
these follow the same trend defined from the Hipparcos stars 
of Section [2] 

We note that the discontinuity is as clear in the purely 
observational plane most sensitive to temperature and lumi- 
nosity (Mk s — (V — K)) as it is in the Mb i — T c g plane: we 
thus consider the temperature-luminosity discontinuity to 
be real. It occurs at 4200 — 4300 K, appearing as a plateau 
in the temperature-luminosity plane, with the luminosity of 
the M dwarfs holding fairly steady even as their tempera- 
tures decrease. For this to occur, the radii of the M dwarfs 
must be increasing again, rather than falling monotonically 
going down the main sequence. 

In Figure 1121 we compare our stars with the very low 
mass star evolutionary models of Baraffe et al. (1998). We 
have adopted 5 Gyr old isochrones, although the evolution 
of the lower main sequence is practically insensitive to the 
age for Mb i > 5.4 (e.g. Casagrande et al. 2007). Figure 
1121 shows that these models do not appear to reproduce the 
discontinuity in the main sequence at 4200 — 4300 K. 

For further insight in the problem, in Figure [13] the 
observed radius-luminosity, radius-mass and the mass- 
luminosity relations are compared with the theoretical pre- 
diction from the same Baraffe et al. (1998) models. We have 
estimated the masses of our stars using the empirical K 
band mass-luminosity calibration of Delfosse et al. (2000) 
which applies for Mk 4.5. For brighter luminosities we 
have used the empirical relation in K band from Henry & 



McCarthy (1993). The use of two different calibrations — 
which are however fully consistent — might be responsible 
for same small offset, but the overall trend is well defined. 
We also convert our 2MASS photometry into the CIT sys- 
tem (Carpenter 2001) before applying the aforementioned 
empirical calibrations. In panels (a) and (b) of Figure [T3] it 
is obvious that current models underestimate the radii of the 
M dwarfs by 15 — 20 percent, as already noticed by several 
other authors (see Ribas 2006 for a review). Such a definite 
conclusion has been obtained using double-lined eclipsing 
binaries. Our study reinforces the finding and confirms its 
existence also to single field stars (Berger et al. 2006). For 
the mass-luminosity relation in Figure [13] (c) , the disagree- 
ment between data and theoretical models is less dramatic. 
In particular, going to masses below 0.5 Mq the agreement 
improve considerably as already noticed by other authors 
(e.g. Delfosse et al. 2000). Since such a good agreement be- 
tween the data and the models is not present in the other 
panels of Figure [13] neither in the temperature-luminosity 
plane of Figure 1121 it argues in favour of a scenario in which 
the stars have larger radii and cooler effective temperatures 
than predicted by models, but just in right proportion to 
barely affect the luminosities. 

In what follows we briefly discuss possible mechanisms 
responsible for the radius increase which marks the transi- 
tion from K to M dwarfs. An interesting discussion on the 
disagreement between the predicted and measured radii of 
very low mass stars from eclipsing binaries and interferom- 
etry can be found e.g. in Lopez-Morales (2007). We also 
mention that another discontinuity at cooler effective tem- 
peratures (V — Ic ~ 2.7), i.e. when the M dwarfs become 
fully convective, is already known in literature (e.g. Hawley 
et al. 1996; Clemens et al. 1998; Koen et al. 2002) and we 
do not discuss it here. 

8.1 Mixing- length 

Very low mass stars are a very interesting place to test the 
input physics in stellar models, since below ~ 0.4 Mq (de- 
pending on the metallicity and the inclusion of magnetic 
fields in the models) stellar interiors are expected to become 
fully convective. Their evolution is thus practically insensi- 
tive to the mixing length parameters omlt and the models 
thus are not subject to any adjustable parameter other than 
the helium abundance (which is expected to be solar scaled) . 
For this reason, very low mass models do not need to be cal- 
ibrated on the Sun. 

The Baraffe et al. (1998) models are computed assuming 
a mixing length qmlt = 1, quite different to values of 1.5 — 
2 which are typically adopted for the Sun, and it is this 
which leads to the difference between the models and data 
for the G and K dwarfs in Figure 1121 In fact, if a solar 
calibrated model is used (Baraffe et al. 1998, continuous 
line), the agreement for those stars becomes excellent. We 
have already extensively tested theoretical models for G and 
K dwarfs in our previous paper (Casagrande et al. 2007) and 
so we focus here on the M dwarfs. 

A possible solution to the discontinuity in the HR di- 
agram could be a rapid decrease of the mixing-length as a 
function of stellar mass, although this would keep rather 
unaffected the lower part of the HR diagram in Figure 1121 
where theoretical isochrones would still remain offset with 
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Figure 12. HR diagram in different planes. Only stars with Hipparcos parallaxes better than 15 percent are used. Overplotted are also 
the theoretical isochrones from Baraffe et al. (1998) for [M/H] = 0.0 (dotted) and [M/H] = —0.5 (dashed) for «mlt = 1- The continuous 
line is a solar calibrated model with »mlt = 1-9- For all the isochrones the age used is 5 Gyr. In the last panel both our sample of stars 
and the theoretical isochrones have been plotted adopting Mg a [Q = 4.74. Notice that the metallicities for the stars plotted here are 
either from high-resolutions spectroscopy (Casagrande et al. 2006) or from the Bonfils et al. (2005) calibration. 



respect to the observed stars. Since the mixing-length de- 
scribes the efficiency of the convection, any physical mecha- 
nism inhibiting convection (like magnetic activity discussed 
in Section f8 . 2f) can be phenomenologically mimicked by de- 
creasing the mixing-length (Chabrier, Gallardo & Baraffe 
2007). Very interestingly, there are indications of a possi- 
ble dependence of the mixing-length with mass from mod- 
elling the components of binaries (e.g. Lebreton, Fernandes 
& Lejeune 2001; Yildiz et al. 2006). Just how viable this 
solution is we regard as an open question. 

8.2 Magnetic activity 

The discrepancy between the predicted and observed radii 
in M dwarfs could be related to the activity level of the stars 
(e.g. Torres & Ribas 2002; Lopez- Morales & Ribas 2005). ft 
is known that active and inactive M dwarfs define two dif- 
ferent sequences in the luminosity-colour (Stauffer & Hart- 



mann 1986) and temperature-radius (Mullan & MacDonald 
2001) plane. Strong magnetic fields are expected to inhibit 
convection, thus giving larger radii for a given T G ff or lower 
T e g for a given radius (Mullan & MacDonald 2001). Alterna- 
tively, it has also been suggested that the larger radii could 
simply be an effect of flux conservation in a magnetic spot- 
covered stellar surface (Lopez- Morales & Ribas 2005). 

The stars with Hipparcos parallaxes plotted in Figure 
1121 are expected to have a very low activity level because of 
the strict selection criteria used in Section [2] This however 
does not exclude the possibility of a large and homogeneous 
spot-coverage since that would not necessarily produce any 
strong variability. To gauge further insight into the problem, 
we have checked that when the stars labeled as variable in 
Koen et al. (2002) are included in Figure [12] they define the 
same trend shown by the stars of Section [2] Notice though 
that the M dwarfs in Koen et al. (2002) were selected among 
the less variable in Hipparcos. Dedicated studies of active M 
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dwarfs should still be done before reaching a more firm con- 
clusion: at present, although magnetic activity can undoubt- 
edly affect stellar radii, we regard it as unlikely as being 
responsible for the main sequence discontinuity observed at 
4200 — 4300 K. Concede the effect of magnetic field certainly 
becomes more important descending the main sequence and 
its inclusion is likely to be a relevant ingredient also for a 
proper modelling of non-active M dwarfs. 

8.3 Opacity 

The models clearly have difficulty in reproducing the strong 
transition between late K and early M type dwarfs, but 
the disagreement becomes even more pronounced as one 
descends to the bottom of the main sequence. Similar dis- 
agreement was already noticed by Baraffe et al. (1998) when 
comparing their models with a sample of field stars, open 
and globular clusters with metallicities similar to those cov- 
ered in the present study. Since the disagreement is much 
less pronounced in the comparison with metal poor globu- 
lar cluster M dwarfs (Baraffe et al. 1997), the disagreement 
at high metallicity might be ascribed to missing opacity of 
some sort in the models. 

Figures [12] and [13] also indicate that the disagreement 
is more marked for metal rich stars, again suggesting that 
missing opacity sources could be a viable solution (Berger et 
al. 2006) . It is very interesting that the discontinuity occurs 
at a temperature when molecular formation (H2O and TiO) 
starts to be important, again suggesting that opacity is a 
possible culprit. 

8.4 Three characters in search of an author 

We have briefly discussed three possible causes for the 
luminosity-temperature discontinuity in the main sequence 
going from K to M dwarfs. The data in Figure[T2]and[T3]sug- 
gest that the problem is more likely to be related to molec- 
ular opacity, than being structural, but considerable further 
work is needed to test those ideas. Simple steps forward to 
confirm or rule out possible explanations would be the anal- 
ysis of late K and early M dwarfs' spectra as well as to 
run the MOITE for a large sample of (magnetically) active 
M dwarfs, to help in searching for correlation between the 
radius discrepancy and activity level or other physical pa- 
rameters (e.g. Lopez-Morales 2007) . Of course, there may be 
no a unique culprit for the radius discontinuity in M dwarfs, 
and only advances in modelling both the structure and the 
atmosphere of these stars will get things right. 



9 CONCLUSIONS 

We have determined the temperature scale of M dwarfs, us- 
ing stars with very accurate multi-band photometry from 
optical to near-IR and the MOITE, a new method which 
exploits the flux ratio in different bands as a sensitive 
indicator of both effective temperatures and metallicities. 
Our proposed temperature scale extends down to T c g ~ 
2100 - 2200 K i.e. to the L dwarf limit (e.g. Leggett et al. 
2002) and above ~ 3000 K is supported from interferometric 
angular diameters. Our metallicities, which are ultimately 
calibrated on Bonfils et al.'s (2005) metallicity scale, are also 



found to be in very good agreement with the latest measure- 
ments from Woolf & Wallerstein (2005, 2006) and Bean et al. 
(2006a, 2006b), even if significant differences in the various 
effective temperature scales still exist. Cool M dwarfs with 
metallicities based on (hotter) Hipparcos common proper- 
motion companions also suggest our metallicities are reli- 
able even below 3000 K, although further data are needed. 
Accurate multi-band photometry for the coolest Hipparcos 
common proper-motion pairs would permit one to firmly ex- 
tend the MOITE to the bottom of the main sequence, thus 
opening this elusive area also to galactic chemical evolution 
investigations. Exoplanets are found around M dwarfs, and 
a uniform metallicity scale for their host stars will also be 
very useful. 

The high quality of our data allows us to identify a strik- 
ing feature which marks the transition from K to M dwarfs, 
which appears to be due to an increase in the radii of the 
early M dwarfs relative to late K dwarfs. We have compared 
our sample of stars with theoretical isochrones for low mass 
stars and find that such a feature is not predicted by the 
models, substantially confirming the disagreement already 
noticed in the case of eclipsing binaries. Possible explana- 
tions including the effect of magnetic fields and molecular 
opacity have been discussed. 

This work also highlight the potentiality of high ac- 
curacy multi-band photometry in determining fundamental 
stellar parameters and identifying fine details in the HR di- 
agram. In particular, the MOITE will hugely benefit from 
the existing infrared (2MASS, DENIS) and forthcoming op- 
tical surveys like SkyMapper (Keller et al. 2007), Pan-Starrs 
(Kaiser et al. 2002) and LSST (Claver et al. 2004) which will 
provide accurate and homogeneous multi-colour and multi- 
epoch photometry for a large number of stars. 
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APPENDIX A: MOITE, TECHNICAL DETAILS 

We use the Phoenix grid of synthetic spectra presented in 
Section[3]to bootstrap the MOITE. We assume log(g) = 5.0 
throughout but we have tested that a change of ±0.5 dex 
in the assumed surface gravity has negligible effect on the 
results. 

For any given star in our sample, we first use the 
T c ff : (R — I)c calibration of Bessell (199f) to obtain an 
initial estimate of the effective temperature T e s t o- We then 
interpolate over our grid of Phoenix model atmosphere to 
compute the flux missing from our multi-band photometry 
and reconstruct the bolometric flux on the Earth. At each 
n-iteration a new T e g >n is obtained — according to equation 
(0 or © — until |T e fr,n - Teff,n-i| < 1 K and the final solu- 
tion is thus found. The rationale that motivates the choice 
between equation (flj or @ will be discussed in the following 
of this Appendix. Notice that at each iteration the estimate 
of the monochromatic and bolometric fluxes also improve 
because of the improved effective temperature used to in- 
terpolate over the grid of model atmosphere. The quantity 
Robs thus is not exactly constant, but it depends — quite 
weakly, indeed — on the improved estimate of the effective 
temperature obtained at each step. 

To interpolate over the grid of model atmospheres, both 
T c ff,n and [M/H] are needed. This is only possible for our 
118 M dwarfs with metallicities obtained from the Bonfils et 
al. (2005) calibration (Section 12. 3p . For the remaining stars 
[M/H] is estimated with the technique presented in Section 
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The behaviour of Rtheo in different bands for various 
metallicities and effective temperatures is shown in Figure 
21 In the infrared, Rtheo increases monotonically with in- 
creasing T e ff above ~ 4000 K. If R oos is greater (smaller) 
than Rtheo, at each iteration T e ff,n increases (decreases) un- 
til it converges to its limiting value. In fact, let us consider 
the following case: 
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and rearranging to highlight the result 
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The case Robs < Rtheo can be similarly proven to give 

Tcff.n < T ff jn _l. 

The technique thus converges quickly above T c g about 
4000 K, even with quite poor initial estimates of the effec- 
tive temperature, as it can be more readily understood by 
looking at the sketch of Figure IA1I Going to cooler effec- 
tive temperatures, a given Robs intersects Rtheo twice, i.e. 
there are two possible iterative solutions. In the case of a 
very cool star (say below 3500 K in the example of Figure 
lAip . if -Robs is greater than Rtheo, at each iteration T c ff, n 
always increases to the solution with highest effective tem- 
perature. Similarly, if R bs is smaller than Rtheo then at each 
iteration T c n in continues decreasing without reaching a solu- 
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Figure Al. Schematic representation of the degeneracy in the 
temperature solution in a given band using the IRFM for cool 
stars. Arbitrary units are used on the j/-axis. At low temperature, 
Robs intersects Rtheo twice, producing two solutions: SI and S2 
(filled circles). If the IRFM starts from A, then R a bs < Rtheo an< I 
at each step the new temperature estimate decreases, and diverges 
away from the desired solution, SI. In case B, R a b s > Rtheo 
and at each step the new temperature estimate increases until it 
reaches the hotter solution S2 and terminates. In neither case is 
the cooler solution SI found. In case C, R bs < Rtheo an d the 
IRFM also converges onto solution S2. The use of equation JTJ in 
the IRFM thus finds only the hotter of two temperature solutions 
and can not be used for very cool stars. A technique which finds 
the appropriate solution over many bands has been developed in 
this paper, and the technical details are discussed here. 



tion. Thus, equation ([T]) cannot be used for the cooler stars 
since it only finds one temperature (the hotter one) or none 
at all. An alternative approach to overcome this limitation 
would be to sample the entire Rtheo space, find the two effec- 
tive temperatures that minimize \Rtheo — Robs\ and choose 
the proper solution. Other then being more computationally 
demanding, at cool temperatures the two minima are quite 
shallow in the infrared since Rtheo inverses smoothly and it 
might not be obvious which one of the two solutions must be 
chosen. On the contrary, the use of the flux products when 
Rtheo increases with decreasing T c ff allows to converge at 
cool temperatures via equation (|3j)- 

It is clear that in any given band, when Rtheo increases 
with increasing effective temperature equation |T} must be 
used, whereas when Rtheo increases with decreasing effective 
temperature equation ([3]) must be used. Only when Rtheo 
inverts any dependence on the effective temperature is lost. 
Fortunately, using multi-band photometry, at any given T e s 
is always possible to find one or more bands for which Rtheo 
has a well defined behaviour, i.e. is either monotonically in- 
creasing or decreasing, discarding bands for which the de- 
pendence is practically flat. 

Above 4000 K we use equation |T} to estimate T c ff. n 
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from J, H, Ks photometry identically to Casagrande et al. 
(2006). Note that the choice between equation (JTJ) or (|3j> is 
important to correctly converge to T e fr in - However, when the 
final solution is found, equation ^ or © returns effective 
temperatures that agree within ~ 5 K in any given band. For 
this reason, once a solution is found it is possible to have an 
estimate of the effective temperature from all the other £ 
bands (T"f). The use of both optical and infrared colours 
is crucial for estimating the metallicities, as we explain in 
Section 14.31 In the optical we estimate from V,Rc,Ic 
photometry; we did not use the U and B bands since these 
colours are not available for all the stars and their theoretical 
modeling is also more uncertain. 

For 2500 < T cff < 4000 K we use equation (JTJ| to esti- 
mate Teff, n from J band and equation to estimate T e w,u 
from V and Rc band. Below 2500 K we also use equation 
((3| to estimate T e a in from Ic- We then average the results 
obtained in these bands for the next iteration. Again, when 
a solution is found we compute the effective temperatures 
Tj predicted by all the colours (with the exception of U and 
B bands, as we already said). For the most metal poor stars, 
Ic and J bands flatten out at very cool temperatures. As 
suggested by Figure [4] we have also implemented a more re- 
fined approach, to ensure that we always use as many bands 
as possible for which Rtheo is expected to monotonically in- 
crease or decrease. Our code, written in IDL, is available 
upon request. 

It is important to note that our technique assigns equal 
weight to each of the bands in the converge to an effective 
temperature. It might plausibly be improved by identifying 
the bands which are more sensitive to effective temperature 
and those to metallicity in determining these parameters. 
In this sense, the IRFM can be regarded as a more elegant 
technique to determine T c s, since it works in the Rayleigh- 
Jeans part of a spectrum and is not much affected by the 
metallicity. However, below ~ 4000 K practically all bands 
start to show considerable dependence on the metallicity, 
the only exception being J band (Figure [3J. The IRFM is 
not quite metallicity independent in any case, since the re- 
construction of the bolometric flux from multi-band pho- 
tometry still depends on the metallicity used to interpolate 
in the grid of model atmospheres. We have looked at two 
scenarios in which MOITE may need to perform, such as 
only IR or optical photometry being available. 

Firstly, we have checked whether any major difference 
arises by using only J band to determine effective temper- 
atures. The behaviour of the flux ratio in such band is in 
fact expected to be quite unique, with very little metallicity 
dependence and always increasing as function of T e ff (Figure 
[4|. We did not find any considerable improvement, but only 
a mean temperature difference of 48 ± 57 K, which we think 
it stems from the zero-point uncertainties in the J band 
absolute calibration (uncertainty which instead average out 
using many bands). In addition, relying on one band means 
the technique is much more exposed to the quality of the 
photometry in that band. 

Secondly, we have studied how the convergence in T e g 
is affected using only BVRc colours. According to Figure [4] 
in the temperature range expected for our stars, it should 
be possible to use optical colours only. We regard the metal- 
licity as a fixed known parameter and we compute the tem- 
perature difference with respect to T e g obtained using both 



optical and infrared colours. We have tested also the dif- 
ference when infrared colours are still used to recover the 
bolometric flux but not for converging in T e e and when the 
infrared colours are not used at all. Summarizing, the mean 
difference is of order 15 ±35 K. This result is very reassuring 
and also makes the technique promising to be used for M 
dwarfs for which only optical colours are available. 

Therefore, at present, the use of all or only of some 
optical and infrared bands seems to return reliable and con- 
sistent results. We plan to further test our findings in forth- 
coming studies, by addressing specifically the sensitivity of 
different spectral bands to effective temperature and metal- 
licity and eventually refine the technique presented here. 



APPENDIX B: THE VEGA AND SIRIUS 
ABSOLUTE CALIBRATION 

In this work we have updated the absolute calibration of 
Vega in the optical by adopting the new reference spectra of 
Bohlin (2007) rather than that of Bohlin & Gilliland (2004) 
and which is expected to be accurate within ~ 1 percent 
in the range 3200 - 10000 A. In the infrared the absolute 
calibration of Vega is kept the same as in Casagrande et al. 
(2006), which is based on Cohen et al. (2003). In terms of 
zero-points, the updated fluxes of Bohlin (2007) corresponds 
to changes of few millimag and affect the derived Teff by 
10 K, thus confirming the results obtained in Casagrande et 
al. (2006). 

For some of our stars, we have also U photometry (Sec- 
tion [2]T|. For Vega we adopt U = 0.02 and the same mag- 
nitudes as in Casagrande et al. (2006) for the other bands 
(i.e. BV {RI)c J H K s) ■ U filters has proven rather difficult 
to standardize (e.g. Bessell 1986, 1990b) and adjustment of 
the U zero-points for different temperature ranges has also 
been discussed (Bessell, Castelli & Plez 1998). U,U — B and 
B — V colours are computed according to the prescription 
in Bessell (1990b). 

The Vega zero-points and absolute calibration thus 
seem now firmly established in the optical (Bohlin & 
Gilliland 2004; Bohlin 2007), but some uncertainties (that 
however do not exceed few percents) still remain in the in- 
frared due to its pole-on and rapidly rotating nature. The 
IRFM and MOITE temperature scales are intimately related 
to the adopted infrared zero-points and absolute calibration. 
The possibility of basing our technique on a different photo- 
metric system and standard star is a valuable sanity check 
to the proposed temperature scale. 

The SAAO JHK photometric system was established 
by Glass (1974) and its accuracy and zero-points refined 
and improved over the years by Carter (1990) and Carter & 
Meadows (1995). Since Vega is unobservable in the Southern 
hemisphere, the zero-points of the SAAO JHK photometric 
system are based on 25 main sequence stars ranging from 
spectral type Bl to A7 (Carter, 1990). Sirius is often chosen 
as an complementary or alternative standard to Vega (e.g. 
Cohen et al. 1992). Its observed magnitudes and colours in 
the SAAO JHK photometric system are given in Table |BT1 

Since no absolute flux measurements are available for 
Sirius, Cohen et al. (1992) decided to absolutely calibrate 
a Kurucz (1991) Sirius model with respect to Vega, by us- 
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Table Bl. Observed magnitudes for Sirius in the SAAO JHK 
system. 



J 


H 


K 


Ref. 


-1.387 


-1.378 


-1.369 


Bessell et al. 1998 



ing observed magnitude difference between Vega and Sirius 
in different near and mid infrared bands. Their resulting 
angular diameter for Sirius was 9 = 6.04 mas, 0.9ct larger 
than the direct measurement (corrected for limb-darkening) 
9 = 5.89 ± 0.16 mas by Hanbury Brown, Davis & Allen 
(1974). Recently, new interferometric measurements have 
become available for Sirius. Davis & Tango (1986) obtained 
9 = 5.93 ± 0.08 (when updated limb-darkening coefficients 
are used, see Kervella et al. 2003), while Mozurkewich et al. 
(2003) found 6 = 5.993 ± 0.108. All these direct measure- 
ments however were obtained at optical wavelength, where 
the limb-darkening corrections are larger and more difficult 
to assess. Recently Kervella et al. (2003) have observed Sir- 
ius in the near infrared, where the limb-darkening correc- 
tions are much smaller, obtaining 9 = 6.039 ± 0.019 mas, in 
superb agreement with spectrophotometric value of Cohen 
et al. (1992). 

We absolutely calibrate Sirius by scaling its latest Ku- 
rucz (2003) synthetic spectrum with the angular diameter 
measurement of Kervella et al. (2003) and therefore indepen- 
dently of any consideration about Vega. The corresponding 
effective wavelength and absolute calibration in the SAAO 
JHK filters are reported in Table lB2l The error in the angu- 
lar diameter given by Kervella et al. (2003) implies an uncer- 
tainty of only 0.6 percent in monochromatic absolute fluxes. 
We adopt a more conservative approach, by taking the stan- 
dard deviation from all the aforementioned interferometric 
measurements: these give an uncertainty of 0.066 mas that 
translates into an uncertainty of circa 2 percent in fluxes, in 
good agreement with the global uncertainty of 1.46 percent 
estimated by Cohen et al. (1992). Also, the fact the domi- 
nant H opacity in A stars is expected to be well understood 
gives confidence on the adoption of a synthetic spectrum. 

We have run the MOITE for stars with SAAO JHK 
photometry ("Section 12. 2p . adopting the zero-points and ab- 
solute calibration of Vega in UBV(RI)c and of Sirius in 
JHK fTable[BT1and[B2l. The difference with respect to the 
use of the 2MASS JHKs photometry and absolute calibra- 
tion (Cohen et al. 2003) is negligible, thus confirming the 
adequacy of the absolute calibration adopted in this work 
and in Casagrande et al. (2006) . The mean difference in T e g 
is 9 ± 3K (a — 25 K), and in both bolometric luminosity 
and angular diameters is well below 1 percent. The results 
provided by the adoption of the absolute calibration of Vega 
or Sirius are therefore identical within the errors. 
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Table B2. Absolute calibration and effective wavelength of the ground-based SAAO JHK photometry of Sirius. Quantities tabulated 
correspond to the definition of the zero magnitude in each filter. 



Band 


^eff 


Monochromatic Absolute Flux 


Uncertainty 




A 


erg cm~ 2 s _1 A _1 


erg cm' 2 s _1 A _1 


J 


12044 


1.176e-09 


2.570c-ll 


H 


16282 


4.079e-10 


8.916e-12 


K 


22004 


1.367c-10 


2.988e-12 



The Kurucz model adopted for Sirius has T eff = 9850 K, log(g) = 4.3, [M/H] = +0.4 and microturbolent velocity £ = km s _1 . The 
same formalism adopted in Casagrande ct al. (2006) is used. Notice that the SAAO JHK photometer is equipped with a InSb detector 
and therefore in generating fluxes from model atmosphere energy-integration is the most appropriate. 



